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ABSTRACT 

We present detailed estimates of "type-I" migration rates for low-mass proto-planets 
embedded in steady-state T-Tauri a-disks, based on Lindblad torque calculations 
ignoring feedback on the disk. Differences in migration rates for several plausible 
background disk models are explored and we contrast results obtained using the 
standard two dimensional formalism of spiral density wave theory with those obtained 
from a simple treatment of three-dimensional effects. Opacity transitions in the 
disk result in sudden radial variations of the migration rates. Regions with minimal 
migration rates may be preferred sites of gravitational interactions between proto- 
planets. Three-dimensional torques are significantly weaker than two-dimensional ones 
and they are sensitive to the surface density profile of the background disk. We find 
that migration times in excess of runaway envelope accretion times or T-Tauri disk 
lifetimes are possible for Earth-mass proto-planets in some background disk models, 
even at sub-AU distances. We conclude that an understanding of the background disk 
structure and "viscosity" , as well as a proper treatment of three-dimensional effects in 
torque calculations, are necessary to obtain reliable estimates of "type-I" migration 
rates. 

Subject headings: accretion, accretion disks — planetary systems: 
proto-planetary disks — planetary systems: formation — waves 



1. Introduction 

Since the first discovery of a planet orbiting a nearby Sun-like star about a decade ago (Mayor 
& Queloz 1995), radial velocity surveys of a few thousands of our closest stellar neighbors have 
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uncovered more than a hundred such planets (Marcy & Butler 1998; Gumming, Marcy &; Butler 

1999; Marcy, Cochran &: Mayor 2000). Transit searches, which have confirmed the gaseous giant 
nature of tlicsc planetary companions (Charbonneau et al. 2000), promise to uncover a large 
number of additional extrasolar planets in the future (Udalski ct al. 2002; Konacki et al. 2003). 

In addition to the large eccentricities of many of the planets discovered to date, one of the 
most striking features of this new population is the existence of a subset of extrasolar planets 
orbiting very close to their parent star (see, e.g., the extrasolar planet almanac^ and encyclopedia"^ 
for a census). It is generally accepted that these close-in extrasolar giant planets could not have 
formed in situ, but instead must have migrated from large distances, where conditions for planet 
formation are more favorable (see, e.g., Lin, Bodenheimer & Richardson 1996). 

There are two physically distinct scenarios considered for the formation of gaseous giant 
planets (see, e.g., Ruden 1999; Wuchterl, Guillot & Lissauer 2000 for reviews). In the standard core 
accretion scenario, dust grain accumulation in a proto-planetary disk proceeds to the formation 
of planetesimals, which later assemble to form proto-planetary cores. Above a mass threshold of 
a few tens M®, runaway accretion of the surrounding gas onto the proto-planetary cores occurs, 
leading ultimately to the formation of giant planets made of gas for the most part (Safronov 1969; 
Mizuno 1980; Pollack et al. 1996; see also Rafikov 2003). Cores of giant planets form more easily 
beyond the "snow line," whose location is typically at AU distances from the parent star (Hayashi 
1991; Boss 1996; Sasselov & Lccar 2000). In the alternative disk instability scenario, a rather 
massive proto-planetary disk is subject to gravitational instabilities leading to the direct formation 
of (possibly multiple) Jupiter-sized objects (Cameron 1978; Boss 1997; 2000, Armitage & Hansen 
1999; Mayer et al. 2002; Lufkin et al. 2003). In the present study, our interests are more focused 
on the core accretion scenario because the evolution of embedded proto-planets is of general 
interest, whether or not gravitational instabilities are important for giant planet formation. 

In both planet formation scenarios, inward migration is required to explain the existence of 
the population of close-in extrasolar giant planets. Here again, two physically distinct scenarios 
for migration have been put forward. Spiral density waves launched by a proto-planct embedded 
in a gaseous disk apply a negative torque on the proto-planet, which migrates inward as a result of 
orbital angular momentum losses (Goldreich &; Tremaine 1980, henceforth GT80; Lin & Papaloizou 
1986; Korycansky & Pollack 1993, henceforth KP93; Artymowicz 1993a,b; Ward 1997a). On 
the other hand, later in the system's evolution, it is also possible for a planet embedded in a 
sufficiently massive planetesimal disk to experience inward migration via repeated planetesimal 
scattering events (Murray et al. 1998). Our interest in the present study is in the gaseous 
migration mechanism, partly because it surely affects the precursors of gas giants, at least within 
the core-instability scenario, whether subsequent planetesimal scattering migration occurs or not. 
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It has long been recognized that the standard core accretion and gaseous migration 
scenarios face a serious difficulty when combined. While estimated migration times for massive 
proto-planctary cores located at ~ AU distances from their parent star are ~ 10^ years (see, e.g., 
Ward 1997a; KP93), estimates for the time required to accrete a gaseous envelope are typically in 
the range ~ 10^-10"^ years (see, e.g., Pollack et al. 1996). Short inward migration times imply that 
cores would be accreted by the central star before they could build up any substantial gaseous 
envelope, and more generally, they pose a serious threat to the survival of any planetary system in 
formation (see, e.g., Ward 1997b). 

Robust calculations for proto-planetary core migration times are notoriously difficult to 
achieve, however. While a forming planet with a mass ^ lO-lOOM® may open a gap and see its 
subsequent migration tied to the slow viscous evolution of the gaseous disk ("type-II" migration), 
a lower-mass proto-planetary core migrates in, relative to the gaseous component ("type-I" 
migration), at a rate which strongly depends on the exact disk conditions in the vicinity of the 
proto-planet. In that respect, most calculations of "type-I" migration rates to date were made 
with rather idealized models for the background disk. Both linear torque calculations (Ward 
1997a; KP93) and fully non-linear hydrodynamical simulations (Kley, D'Angelo & Hcnning 2001; 
Nelson & Bcnz 2003a, b; D'Angelo, Kley &: Hcnning 2003; Lufkin et al. 2003) usually assume power 
law profiles for the disk properties, in models often reproducing the minimum mass solar nebula 
(MMSN) characteristics. This is one of our main motivations to study "type-I" gaseous migration 
with more realistic proto-planetary disk models. While there are significant uncertainties regarding 
the nature of angular momentum transport in these disks (Gammie 1996; Glassgold, Najita & 
Igea 1997; Fromang, Terquem & Balbus 2002; Fleming & Stone 2003; Matsumura & Pudritz 2003) 
or their overall structure (e.g., Chiang & Goldreich 1997), we believe that a study of how these 
uncertainties may affect migration time estimates is useful at this time. 

Another source of uncertainty in calculations of "type-I" migration rates comes from 
limitations of spiral density wave theory itself. The theory has been initially developed for 
infinitely-thin, two-dimensional disks (Goldreich & Tremaine 1980), but three-dimensional effects 
are important (Tanaka, Takeuchi & Ward 2002) . Finite eccentricities and inclinations of embedded 
proto-planets (Papaloizou &; Larwood 2000; Goldreich & Sari 2003; Ogilvie & Lubow 2003) and 
the potentially magnetized (Terquem 2003) and turbulent nature (Winters, Hawlcy & Balbus 
2003; Papaloizou &: Nelson 2003; Nelson & Papaloizou 2003a, b; Papaloizou, Nelson & Snellgrove 
2003) of the proto-planetary disk are also ignored in the original theory, while they could all have 
important consequences. It is probably fair to say that a consensus has yet to be achieved on 
these various aspects of spiral density wave theory. 

Nevertheless, an important motivation for developing a detailed numerical tool to study 
gaseous migration in proto-planetary disks is the intimate relation, which has become increasingly 
obvious, between planetary formation and migration on the one hand, and disk evolution on the 
other hand. Planets form from disk material, migrate by interacting with the disk and at the same 
time influence strongly the disk long-term evolution by raising torques which affect the disk global 
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angular momentum budget and eventually open gaps in an otherwise smooth gaseous distribution. 

Therefore, it is likely that a general planet formation theory capable of explaining the diversity of 
known systems will require following the strongly-coupled evolution of a potentially large ensemble 
of proto-planets together with the disk in which they formed (see, e.g., Lin & Papaloizou 1986 for 
early calculations of this type). 

In this first investigation of gaseous migration in proto-planetary disks, we focus on "type-I" 
migration of low-mass proto-planets, ignoring the feedback from Lindblad torques on the disk 
structure. The radial structure of our disks is not a simple power law but depends upon opacities 
and other details of accretion. We explore how uncertainties in the structure of proto-planetary 
disks may affect "type-I" migration rates and how sensitive these rates are to three-dimensional 
effects. In §2, we describe the characteristics of our T-Tauri disk model, while in §3 the two- 
and three-dimensional formulations we have adopted for Lindblad torques are summarized. 
Corotation torques, and our reasons for neglecting them, are discussed briefly in §3.3. We present 
our results in §4, we discuss possible limitations and extensions of this work in §5 before concluding 
in S6. 



The present study is exclusively focused on steady-state a-disk models (see §4.1 for some 
justification). In the future, however, we intend to study the coupled evolution of planet-disk 
systems with explicit planetary migration, and this obviously requires being able to follow the 
evolution of the disk's properties with time. With this in mind, we have developed a fiilly 
time-dependent model for T-Tauri a-disks. It is largely inspired by the model developed by 
Hameury et al. (1998) to study the time-dependent evolution of disks around compact objects. 
The main assumption in our models is that disks are geometrically thin, so that their radial and 
vertical structures can effectively be treated separately. 



Following Hameury et al. (1998), we solve the following radial equations for the conservation 
of mass, angular momentum, and energy: 



2. 



Disk Models 



2.1. 



Radial Structure and Evolution 




(1) 
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where S is the surface column density, is the gas radial velocity in the disk, j — (GM^r)^/^ 
is the specific angular momentum of material at radius r in the disk, fix = (GM^/r^)^/'^ is the 
Keplerian angular velocity, is the mass of the central star (taken to be 0.5Mq in the present 
work), Tc is the mid-plane temperature, Cp is the heat capacity at constant pressure, 3ft is the 
perfect gas constant, /j, is the mean molecular weight and v is the kinematic viscosity coefficient, 
for which we adopt a standard a-prescription (Shakura &; Sunyaev 1973). Justifications for the use 
of an a -viscosity to describe angular momentum transport in thin accretion disks have been given 
by Balbus &: Papaloizou (1999) in the case of magnetized turbulent transport and by Goodman &; 
Rafikov (2001) in the case of transport by locally-dissipated planetary torques. The gas motion is 
implicitly described as Keplerian in our formulation of the disk equations. Deviations from strict 
Keplerianity due to finite radial pressure gradients are calculated separately, when required for 
the planetary torque calculations, and deviations from Keplerianity due to disk self-gravity are 
ignored. 

The terms S and J in the first two equations represent additional sources/sinks of mass and 
angular momentum, respectively, for the disk (both are put to zero in the present study). Q'^ and 
Q~ are the local heating and cooling rates per unit surface, respectively. The term H accounts 
for the radial energy flux caused by steep radial gradients, and it is estimated in the framework 
of the a parametrization (see Hameury et al. 1998 for details). This general form of the energy 
equation accounts for the possibility that the disk is out of thermal equilibrium (Q+ / Q~) as 
can occasionally happen in time-dependent studies (see, e.g., Menou, Hameury & Stehle 1999 for 
specific examples in dwarf nova disks). For all the steady-state models presented here, however, 
the radial energy equation reduces to a very good approximation to = Q~ . 

Equations (1) and (2) can be combined to obtain the standard viscous diffusion equation for 
a geometrically thin disk (when S = J = 0; see, e.g., Pringle 1981). The above set of equations is 
solved on an adaptive numerical grid, to guarantee that regions of the disk characterized by sharp 
gradients receive a large number of grid points and thus adequate numerical resolution. A number 
N = 400 grid points has been adopted for all the calculations presented here. The equations are 
evolved with a fully implicit numerical scheme and a variable time-step. Further details can be 
found in Hameury et al. (1998). For the present study, steady-state disk models were obtained in 
the limit d/dt ^ (in practice, a large enough time-step). 

This set of partial differential equations requires a total of six boundary conditions. For a disk 
with fixed inner and outer radii which is fed, at its outer edge, at a rate M = Mout, the boundary 
conditions at the inner edge (n = 1) and the outer edge {n = N) of the numerical domain are: 

e < 1 g cm~^, 
Tin (0.01 AU), 

0, 

Mout, 



r\n=l = 

n=l — 

or 

ML=N = 
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r\n=N = rout (200 AU), 
dT 

-Q^\n=N = 0. (4) 

The precise value of S|„=i is unimportant as long as it is small compared to values in the bulk of 
the disk. 



2.2. Vertical Energy Transport 

Contrary to Hameury et al. (1998), we do not explicitly solve for radiative-convectivc energy 
transport in the vertical direction. Instead, we calculate the local cooling rate Q~ = 2aT^Q (for 2 
disk faces) with a one-layer version of the radiative transfer model of Hubeny (1990). According 
to this model, the disk intrinsic effective temperature, T^g, is related to the disk mid-plane 
temperature, Tc, by 



T-i4 ^ (-^c -^irr) /r\ 

3(rext/2 + l/V3 + l/(3r,bs)) 

where the mid-plane optical depth r = kT,/2, Kcxt is the Rosseland-mean extinction opacity 
(absorption+scattcring), Kabs is the Rosscland-mcan absorption opacity and T[j-j- > in the 
presence of external irradiation. Note that in this formulation, T^g accounts only for the disk 
intrinsic emission (i.e. from viscous dissipation only), while the disk also re-emits all the irradiation 
flux absorbed locally (as measured by Tirr). A detailed discussion of this radiative transfer model 
and of how well it does in the two important optically-thick and optically-thin limits can be 
found in Hubeny (1990). While Hubeny's model requires in principle a Planck-mean opacity to 
be used for Tabs (in the optically-thin regime), we only use Rosseland-mean opacities here (see 
§2.3 for a justification). We do not allow for different opacities when considering the disk's own 
radiation and the hotter radiation from the central star. Adopting this simple prescription for 
vertical radiative transfer is an oversimplification of the problem, which is justified a posteriori 
by satisfactory comparisons to previous studies treating the vertical energy transport problem in 
more details. Our one-layer approximation neglects vertical transport by convection, whose role 
is unclear, especially in the context of magnetized accretion disks, which are already subject to a 
dynamical instability leading to turbulent transport, by construction. We also note that D'Alessio 
et al. (1998) have estimated that convective transport should be unimportant in all but the very 
innermost regions of a typical T-Tauri disk. 

Viscous dissipation is prescribed in the usual form, with a heating rate per unit area 

Q+ = ^i^Jinj,, (6) 

where the kinematic viscosity, u = acgh, is expressed in terms of a dimensionless viscosity 
parameter, a {cg is the local adiabatic sound speed and h is the local scale height). 
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Irradiation by the central star is an important source of heating for the disk at large distances. 

Rather than solving self-consistently for the magnitude of irradiation based on a detailed model 
of the disk shape and structure, as in Dubus et al. (1999), wc parametrize it with a simple 
geometrical model. Assuming that the irradiation source is point-like (a reasonable assumption 
far away from the stellar source), the irradiation flux is given by 

where a is the Stcfan-Boltzmann constant, is the stellar luminosity, e is the disk albedo and 
Tirr should go to zero if the last term in parenthesis becomes negative (disk self-shadowing; see, 
e.g., Dubus et al. 1999). For e = 0.5, an aspect ratio h/r ~ 0.1, a geometric factor ~ 0.1 (last 
term in parenthesis in Eq. [7]) and L* ~ Lq, this yields an irradiation temperature 

which is directly used in equation (5). This prescription artificially exposes disk regions 
with dlog /i/dlog r < 1, that would normally be self-shadowed, to irradiation from the central 
star, but without serious consequences as we shall see below. Although an improved treatment 
of disk irradiation would lead to better accuracy, this simple prescription should be sufficient for 
our purposes in view of other large uncertainties inherent to structural models of proto-planetary 
disks. 



2.3. Opacities and Equation of State 

We assume that dust (when present) is fully mixed with the gas. Analytic expressions 
for the Rosseland-mean extinction opacities are taken from Bell et al. (1997; see also Henning 
& Stognienko 1996). Pollack et al. (1994) have found that Planck-mean and Rosscland-mean 
extinction opacities are comparable for mixed dust. This justifies our use of Rosseland-mean 
opacities in Eq. (5) even in the optically-thin regime. Following Pollack et al. (1994; see their 
Fig. 4b), we also adopt /tabs = /^ext) since we have found that this is adequate for the range of 
temperatures at which typical T-Tauri disks are optically-thin. 

The equation of state adopted is that of a perfect gas. The mean molecular weight and 
adiabatic index are typically /x^jj 2± 2.4 and F 2± 7/5 throughout the disk, since hydrogen 
is everywhere predominantly in molecular form. Other thermodynamical quantities are self- 
consistently calculated as a function of temperature and density at the disk mid-plane when 
needed (see Hameury et al. 1998 for details). 
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3. Planetary Torques 

In this study, we focus on the low-mass pcrturber hmit and we assume that the feedback on 
the disk due to the presence of embedded proto-planets can be neglected. We neglect corotation 
torques in our calculations, and assume for simplicity that all the planets are on circular orbits, 
without inclination relative to the disk mid-plane. Linear spiral density wave theory has been 
elaborated by a number of authors (GT80; KP93; Artymowicz 1993a, b; Ward 1997a) and we rely 
heavily on these works in our formulation of the problem. Thus, we calculate the "torque density" 
representing waves excited at Lindblad resonances, which are located at radii Vm where the disk's 
angular velocity Q(r) and epicyclic frequency ^(r) are related to the angular velocity Jlp of the 
planet by m [ft{rm) — ^p] = K,{rm) for azimuthal harmonics m. 



3.1. Torque Density 

Following Ward (1997a), the torque density (torque per unit radius) of waves excited in the 
disk by a planet of mass Mp = ixM^ orbiting a star of mass M* at radius Vp is 



dT ■ 



= sign(r - Tp) . 



,4-2 



(9) 



The surface density, S, and epicyclic frequency, k (not to be confused with opacity), are 
evaluated at the radius of interest, r, and ^ = mc^/rK, where Cg is the adiabatic sound speed. This 
expression incorporates the torque cutoff via the functional dependence of £,{r), ip{r) and m{r) 
(see below). 



Note that 



(o - UpY 



(n-Qp)^ -c2/r2' 



(10) 



The angular velocity and epicyclic frequency in the gaseous disk are slightly different from 
their Keplerian values due to finite radial pressure gradients: 



where Ciso is the isothermal sound speed. We must include the pressure terms in the disk 
rotation curve when we compute the torque density, even if we neglect them in the equations 
for the disk evolution (§2.1). The angular velocity of the satellite, Q,p, has the Keplerian value 
^GM^/r^, and the small difference Qp — fi(rp) is the most important cause of the imbalance 
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between the torques exerted on the interior and exterior parts of the disk, and hence of the planet's 
migration. Because of the extra radial derivative, k is even more sensitive to pressure effects.^ 

The dimensionless azimuthal wavenumber m is treated as a continuous function of radius. It 
is computed from the WKBJ dispersion relation by setting k ^ m/r (since the radial wavenumber 
vanishes at the Lindblad resonances where the waves are excited): 



m(r) 



(Jl - J^p)2 - c2/r2 



1/2 



(12) 



Where equation (12) predicts imaginary m, dT/dr 0. Disk self-gravity is neglected in 
equation (12). 



The dimensionless satellite forcing function is 



TT 


1 


""l/2 


2 


m 


dp 



(13) 



where P = r/rp. The factor Q/k docs not appear in the original formulation of Ward (1997a) 
but it is needed so that equation (13) reduces to equation (18) of Goldreich & Tremaine (1980), in 
the limit ^ ^ 1. The standard definition of the Laplace coefficient is 



bTM = -[ . ■ (14) 

^/^^^ ttJ ^/l- 2/3 cos + ^2 ^ ^ 



This is exact for an infinitely-thin two-dimensional disk, as studied, for instance, by Ward 
(1997a), and this is the formulation we adopt for 2D Lindblad torques. 

To account for three-dimensional effects, we replace the above expression with the following 
approximation, which includes a correction for the thickness, h, of the disk 

bT/^m ^ ^ ^0 (m^/3-2 + /3-i +hyrrp^ , h^^, (15) 

where Kq is the modified Bessel function of the second kind of order 0: Kq{z) « — Inz for 
< z <^ 1 and Kq{z) ^ e~^y/TT/2z for z S> 1. Following Goldreich & Tremaine (1980), we have 
made use of the integral representation 

oo 

f cos{zt) 



. , , f coslzt) , 



^There is some danger that may be formally negative near the edges of gaps, which are not considered in the 
present study. 
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(Abramowitz &; Stegun §9.6.21). The appearance of a term proportional to in equation (15) 
comes from expressing the gravitational potential of the perturbing planet as 

= - , . (17) 



while h is set exactly to zero in the two-dimensional theory. 

With this expression, we account directly in the perturbing potential for the fact that the disk 
material subject to the gravitational influence of the proto-planct is actually distributed vertically 
over a length scale typically ~ h, the disk scale height. It is clear that this formulation accounts 
only approximately for three-dimensional effects. In a detailed study, Tanaka et al. (2002) find 
that the exact vertical distribution of mass in the disk enters explicitly a three-dimensional 
formulation of spiral density wave theory (see their Eq. [28], which is derived under a specific 
vertically isothermal assumption). As we emphasized earlier, the structure of proto-planetary 
disks, including their vertical structure, is not well known (e.g., Fromang et al. 2002; Fleming & 
Stone 2003; see also Chiang &; Goldreich 1997). In this context, it is unlikely that three-dimensional 
torques can be calculated very accurately. This is our main motivation for adopting a simple 
formulation of three-dimensional effects, that should nonetheless capture their most salient 
features. We note that Papaloizou & Larwood (2000) have adopted a similar formulation, with a 
vertical "softening length" varying from 40% to 100% of a disk scale-height, h (see also Ward 1989 
for a different treatment of vertical thickness). It is unclear what value of the softening length is 
most appropriate or best reproduces the results of Tanaka et al. (2002) in general, but most cases 
of interest should be covered by a treatment intermediate between our two- and three-dimensional 
prescriptions. 

The expression for tp includes a derivative with respect to (3 = r/vp, and we include the 
implicit dependence on r via h?{r)/r when evaluating this derivative. However, m is treated as a 
constant during the differentiation. To evaluate the term db^i^il^i) / d(3 in equation (13), note that 
dKQ{z)/dz = —Ki{z) (Abramowitz & Stegun §9.6.27), where Ki is the modified Bessel function 
of the second kind of order 1. Therefore 



Since h = ciso/^^, 




(18) 



1 — = b-^^^- + - — • 

dr r \dr rjr 

When /i <C |?' — 'Tpl ^ r^, then ^ <C 1 and reduces approximately to a constant: 

2fi 

~ K^{k/\2A\) + —Ko{k/\2A\), 

K 
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where 2A = dQ/dlnr = (k^ — 40^)/2r2 is the shear rate. Unfortunately, the torque is dominated 
by the region |r — rp\ ~ h, where this approximation cannot be used. In fact, as |r — rp\ 2h/3, m 
diverges and ifj ^ exponentially. The torque density vanishes at \r — rp\ < 2h/3: in this region, 
the flow past the planet is subsonic so there is no wave drag. 



3.2. Implementation of Planetary Torques in Disk Models 

The torque density has a sharply peaked profile, cutting off at a distance ~ h from the 
proto-planet's location (see, e.g.. Ward 1997a). An accurate calculation of the torque density 
therefore requires that radial profiles in the background disk model be adequately resolved on 
length-scales smaller than a disk scale height. This is not trivially achieved because the typical 
grid spacing in our disk models is usually not very much smaller than a disk scale height. 

A first method we explored to calculate the torque density was to make use of the adaptive 
grid. In principle, by requiring that the grid spacing be inversely proportional to the local torque 
density, one should reach adequate spatial resolution in the disk model (in the vicinity of a 
proto-planet) for the torque density to be calculated accurately. In practice, we have found that 
this solution is difficult to implement and often results in numerical instability. 

Instead, we have opted for a sub-grid method. Whenever required for a torque density 
calculation, the radial profiles of relevant quantities in the disk are fitted to cubic splines in 
the vicinity of the proto-planet's location. For the torque density calculation itself, cubic spline 
interpolations of the disk profiles provide estimates of disk properties at distances as close to the 
proto-planet as required (whose location does not in general coincide with that of a grid point in 
the disk model). This implementation effectively decouples the accuracy of the torque density 
calculation from the effective numerical resolution of the background disk model. Generalization 
to multiple proto-planets embedded in the disk is straightforward. 



3.3. Corotation Torques 

In addition to the Lindblad torques, there are torques at co-orbital corotation resonances. 
The latter occur very near the semimajor axis of the planet, rc — ~ 0(/i^/r^). Eccentric planets 
have non-co-orbital corotation resonances as well {e.g., GT80). 

The theory of corotation resonances is in a less satisfactory state than that of Lindblad 
resonances. GT80 derived an expression for the corotation torque due to an individual azimuthal 
harmonic. In the co-orbital case and for m S> 1, their expression reduces to 

rS- » -s,nK„Mi - [^)' , m 

where /3 = r/rp and the sign pertains to the torque on the planet. Pressure forces and the 
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disk thickness were neglected; normally these effects should be unimportant sX m r/h. The 

dominant contribTition to the total torque on an embedded planet, however, comes from large 
m. For the Lindblad resonances, GT80 demonstrated that harmonies m > r/h are suppressed 
by pressure even in a two-dimensional disk. This result has been elaborated by later authors as 
summarized above. Analogous "torque cutoffs" were derived for corotation resonances by Ward 
(1989) using asymptotic methods, and by KP93 using numerical integration of the linearized 
equations. The latter concluded that the net torque (summed over rn) from co-orbital corotation 
resonances is comparable to the residual of Lindblad torques — a result that is suggested by simply 
truncating the sum over eq. (19) at m = r/h. Since all harmonics are proportional to the term 
in square brackets above, however, the corotation torque is very sensitive to the local gradient in 
surface density — even in sign — and to the third derivative of the rotation curve. Coincidentally, 
the torque vanishes in a keplerian disk with the canonical profile of the minimum-mass solar 
nebula = Q (xT, oc r^^/'^). Unfortunately, the torque cutoffs derived by Ward (1989) and KP93 
do not agree quantitatively, even though both groups used linear theory. A subset of KP93's 
results were confirmed in the 2D limit of a three-dimensional linear analysis by Tanaka et al. 
(2002), who treated the special (and relatively tractable) case of disks that are both radially and 
vertically isothermal. 

Linear theory itself is problematic. The angular momentum deposited at Lindblad resonances 
is carried away by a wave flux. But the torque at corotation accumulates in a small amount 
of disk material unless spread by viscous diffusion or by differences in the radial drift rates of 
the planet and the disk (Ward 1989; Balmforth & Korycansky 2001; Masset 2001, 2002). Hence 
the corotation resonances are prone to nonlinear saturation. The larger the gradient of K^/2ilS 
(potential vortieity) at the resonance, the stronger is the torque in linear theory, but the more 
sensitive to saturation. In a disk with a smooth profile, saturation is not likely if the planet's mass 
is much less than both of the gap-opening criteria cited in §4.4. 

In summary, we have neglected corotation torques for lack of a simple and reliable 
parametrization of their effects. 

4. Results 

4.1. T-Tauri a-Disk Models 

Although accretion disks in T-Tauri systems are not expected to be in steady state, a 
quasi-steady assumption may be a reasonable first approximation for all but their outermost 
regions. Indeed, assuming the disk is freely spreading, the global disk evolution should be 
determined by comparatively slow viscous processes in the outer regions, while the inner regions 
are able to follow this evolution via a succession of quasi-steady states (see, e.g., Cannizzo, Lee & 
Goodman 1990; Ruden & Pollack 1991). Using this general property, Hartmann et al. (1998) 
have shown that the observational characteristics of pre-main sequence stars can be successfully 
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interpreted with freely-spreading accretion disk models in which the inner regions go through a 

succession of quasi-steady states with an accretion rate that is slowly decreasing with time. They 
infer mass accretion rates in the range 10~'^-10~^Mq yr~^, with typical values ~ 10~®Mq yr~^ for 
10^ year old T-Tauri star systems. 

Based on this, our reference T-Tauri disk model is a steady-state model with a constant 
accretion rate M = 1O~^M0 yr~^ throughout. We adopt = O.SM© for the mass of the central 
star and a default value of the viscosity parameter a = 0.02. Disk irradiation by the central star, 
when included, follows the prescription given in equation (8). Later on, we will vary various disk 
parameters to study how changes in the background disk model influence the planetary migration 
process. We detail here the properties of our reference disk model because it illustrates rather well 
the basic characteristics of all the disk models we have considered. 

Figure la shows the radial profiles, from 0.01 to 200 AU, of the surface density (solid line), 
Rosscland-mcan optical thickness (short-dashed line), aspect ratio (long-dashcd line) and Toomre 
Qt parameter for self-gravity (dotted line) in the reference model. For each quantity, a profile 
with and without disk irradiation is shown. Note that the profiles very near the inner edge of the 
numerical domain are affected by our choice of boundary condition but they become independent 
of it beyond 0.02 AU or so. The disk is characterized by an aspect ratio rising from ^ 10~^ close 
to the central star to ^0.1 at AU distances and beyond. Disk self-gravity is at most weak, even 
for the coldest outer regions, at r = 200 AU. The disk makes a transition from an optically-thick 
regime to an optically-thin regime at a distance of several AUs from the central star. 

Figure lb shows radial profiles of the midplane temperature (solid line) and internal effective 
temperature, arising from internal viscous dissipation only (dashed line), for the same reference 
model with and without irradiation. The transition from an optically thick to an optically 
thin regime is again easily identified at a distance of several AUs from the central star. In 
the irradiated model, the outer, optically-thin regions of the disk are vertically isothermal (the 
midplane temperature being equal to the irradiation temperature). 

In the outer regions of the non-irradiated model, the midplane temperature (lower solid line) 

is significantly in excess of the effective temperature (dashed line) because of the inefficiency of 
the optically-thin gas at radiating the dissipated energy. It is clear from figures la and lb that the 
effect of disk irradiation is relatively minor and that it only is important for regions of the disk 
beyond 1 AU from the central star. Note also that regions of opacity transitions in the disk can 
already be identified in figure 1 as regions where the slopes of many of the profiles shown change. 
As we will see later, these slope changes have potentially important effects on the planetary 
migration process. 

Our results for the steady-state structure of a typical T-Tauri a-disk are in good overall 

agreement with previous comparable studies (see, e.g., D'Alessio et al. 1998; Papaloizou & 
Tcrquem 1999; Fromang et al. 2002). In particular, the comparison to D'Alessio et al. (1998) 
is quite satisfactory given that these authors treated the disk structure in considerably more 
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details than we did. These successful comparisons vahdate a posteriori a number of simphfying 
assumptions that we made in our treatment of the disk structure (see §2). 

This success should not, however, eclipse the fact that there remains considerable theoretical 
uncertainty regarding the structure of proto-planetary disks. Perhaps the most important source of 
uncertainty comes from our ignorance of the nature of angular momentum transport ("viscosity") 
in the bulk of proto-planetary disks (Gammie 1996; Glassgold, Najita Sz Igea 1997; Fromang et al. 
2002; Fleming &; Stone 2003; Matsumura Sz Pudritz 2003). Over the past few years, it has also 
become clear that the role of a super-heated layer of dust is important to explain the observational 
characteristics of T-Tauri disks (Chiang &; Goldreich 1997), an aspect of the problem that has 
simply been ignored in our treatment. 



4.2. Tests of Lindblad Torque Calculations 

Tests of our implementation of the torque density calculations serve several purposes. We 
want to recover previously established results, guarantee numerical convergence and identify the 
main differences between the classical 2D treatment of spiral density wave theory and our simple 
3D treatment. 

Figure 2 shows results for our 2D implementation, in a format such that they can be 
directly compared to similar results obtained by Ward (1997a), for a disk with fixed aspect ratio, 
h/r = 0.07. The inner, outer and differential Lindblad torques (or, equivalently, integrated torque 
densities) are normalized to the value 

T„^.^,rX{^)\^y\ (20) 
where the subscript p refers to the proto-planet's location. 

The background disk is idealized by assuming that the radial profiles of surface density and 
midplane temperature are exact power laws, with indexes 

^^_dlnE^ ^_ dlnT, ^ 
dlnr ' dlnr 

A comparison between our figure 2 and figure 3 of Ward (1997a) shows that we recover his 
results on both the magnitude of Lindblad torques and their dependence on k and I with high 
accuracy. In particular, we confirm that 2D Lindblad torques are essentially independent of the 
background profile of surface density (as measured by k). 

Figure 3 shows very different results for our simple 3D implementation of Lindblad torques. 
The 3D torques are reduced by a factor ~ 3-5 from the 2D values, and their dependence on 
the background disk profiles is affected. More precisely, there is now a clear dependence on the 
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background surface density profile (as measured by k). In fact, the weakness of 3D Lindblad 

torques and their dependence on the background disk properties are such that a sign reversal 
of the differential Lindblad torque appears possible for sufficiently steep surface density profiles 
(large k) and sufficiently shallow midplanc temperature profiles (small I). 

Wc have confirmed that the results shown in figures 2 and 3 are unaffected by changes of 
numerical resolution in the background disk models, in the range N = 100 to 1000. These tests 
show adequate numerical convergence of our torque density calculations based on a sub-grid 
method (see §3.2). 



4.3. Migration Rates 

Given a proto-planet of mass Mp embedded in a disk at a radius r^, the net torque, Jp, 
experienced by the proto-planet is calculated as a function of the disk properties in the vicinity of 
rp. In the present work, we neglect the feedback on the disk structure caused by the proto-planet 's 
influence. This assumption should be reasonably accurate in the low-mass limit. Rather than 
following the actual migration with time of the proto-planet in the gaseous disk, we estimate the 
local inward migration time of the proto-planet as 



J„ M„./GAI,r„ , , 

Tn.ig ~ 3^ ^ '\ (22) 

Jp Jp 

where we have assumed a low-mass proto-planet again in the expression for the Kcplerian 
angular momentum, Jp. By populating the entire disk with a sufficient number of identical mass 
proto-planets, we obtain a radial profile of local migration times. In what follows, we show how 
such migration time profiles are affected by variations in the background disk model. We have 
found that populating the disk with 200 identical proto-planets, logarithmically spaced from 
Tin = 0.01 AU to Tout = 200 AU was enough to obtain accurate migration time profiles in all cases 
of interest. 

Figure 4 shows radial profiles of inward migration time, Tmig, for Earth- mass proto-planets 
embedded in our reference T-Tauri disk model, with M = 1O~^M0 yr^^ and a = 0.02 (see 
Fig. 1 for detailed disk properties). The solid line shows migration times calculated with the 
3D formalism, while the short-dashed line shows the same result according to the standard 2D 
formalism of spiral density wave theory. As expected, migration times arc significantly longer in 
the 3D formulation. The dotted line shows, for the 3D formulation, that the effects of neglecting 
disk irradiation on the estimated migration times are significant only at distances ^ 1 AU from 
the central star, where the disk structure becomes sensitive to irradiation. In addition, the 
long-dashed line shows the profile of migration times obtained with a modified 3D formulation, 
assuming that the gas is vertically extended over only half a local disk scale height for the torque 
calculation {h 0.5h in Eq. [17]). This confirms that migration times are sensitive to the exact 
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vertical distribution of gas in the disk. 

The most striking features of figure 4 are the sudden variations of Tmig with radius, at 
distances ^ 1 AU from the central star. Due to a stronger dependence on the background disk 
properties (in particular on the surface density profile), radial variations of Tmig are much more 
pronounced in the 3D formulation of Lindblad torques than in the 2D one, even if they arc clearly 
present in both cases. It is also apparent that these variations of Tmig occur at specific locations 
in the disk. They are easily identified as the regions where changes of opacity regime occur in the 
disk (see the corresponding slope changes in radial profiles of E, Tc and other quantities in Fig. 1). 
These radial variations of Tmig (or, equivalently, of the differential Lindblad torque) have their 
origin in differences between the inner and outer disk properties when the proto-planet is located 
in the close vicinity of an opacity transition. In these special regions of the disk, with different 
inner and outer slopes for S and Tc, differential Lindblad torques can take values that are quite 
different from more typical values adopted in the rest of the disk. 

This property of differential Lindblad torques is interesting because it suggests that narrow 
regions (e.g., at 0.1 and 0.25 AU in Fig. 4) corresponding to local maxima of the migration times 
may be the preferred sites of mutual interactions between proto-planets. Clearly, this possibility 
and its consequences can only be investigated by actually following the time-dependent migration 
of a population of proto-planets, and accounting for their mutual gravitational interactions as they 
approach stalling regions such as the ones identified in figure 4. 

Another notable feature of the migration times shown in figure 4 is their overall magnitude. 
In agreement with previous studies (see, e.g.. Fig. 14 of Ward 1997a), we find that the typical 
migration time of an Earth-mass proto-planet located at a few AU distances from the central 
star is of order a few 10^ years, according to the 2D formalism of spiral density wave theory. 
As we already mentioned in §1, this relatively short timescale has been recognized as a major 
difficulty for the standard core accretion scenario for giant planet formation because it may not 
allow enough time for the build up of a substantial gaseous envelope. On the other hand, the 
typical migration times obtained with the 3D formulation arc typically one order of magnitude 
longer than the 2D ones, and regions in the vicinity of opacity transitions are characterized by 
migration times approaching 10^ years for an Earth-mass proto-planet. We view this result as an 
indication that a proper treatment of 3D effects in Lindblad torque calculations, together with a 
detailed background disk model, may be essential to the resolution (even if only partial) of this 
long standing difficulty for the core accretion scenario. 

In what follows, we investigate the effect on migration times of variations in the background 
disk model. Since the 3D formalism is arguably the most appropriate (proto-planetary disks do 
have a finite thickness) and it is more sensitive than the 2D one to details in the background disk 
properties, we chose to show results for migration times calculated with the 3D formulation only. 
We also restrict ourselves to disk models that include the effects of irradiation by the central star 
on the disk structure. 
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Figure 5 shows the effect on migration times of varying the viscosity parameter a in 
the disk model. As a is reduced, the disk mass is increased (at a fixed mass accretion rate 
M = IO^^AIq yr^"*^), so is the differential Lindblad torque, so that migration times are reduced. 
A secondary consequence of varying a is to increase the values of the radii at which opacity 
transitions occur in the disk, so that the regions of maximum migration times are shifted outward 
as well. The relative amplitude of radial variations of Tmig remain comparable even as a is varied. 

Figure 6 illustrates the effect of varying the mass accretion rate M in a background disk 
model with the same viscosity parameter a = 0.02 as in our reference model. This time, the overall 
magnitude of migration times is only slightly affected (disks with larger M are somewhat more 
massive), while the profiles are shifted outward as the mass accretion rate is increased. Finally, we 
note that we have checked directly the linear dependence of T^ig with the proto-planet's mass, 
as expected for this migration of "type I" (without disk feedback). Consequently, even for a 
proto-planet of mass Mp = lOM,^, which in some specific core accretion models is considered as 
a typical threshold mass for the onset of runaway envelope accretion, migration times in some 
regions (at sub-AU distances) of our reference T-Tauri disk model are still in excess of 10^ years. 



4.4. Conditions for Gap Formation 

One of the assumptions made in our derivation of migration times is that the feedback 
torque exerted on the disk by a proto-planet is negligibly small. While this assumption should be 
reasonably accurate in the low-mass limit, we can further quantify it by estimating the mass above 
which the feedback torque becomes strong enough for the proto-planet to open a gap. Following 
Rafikov (2002; see also Lin &; Papaloizou 1986 and Ward &; Hourigan 1989), a proto-planet with a 
mass in excess of the "viscous mass" 

/, \5/2 

My ~ (10.34a) i-^j (23) 

generates a torque larger than the disk's own viscous torque and opens a gap. In addition, 
an independent condition for gap formation is that the proto-planet's Hill radius exceed the local 
disk scale-height, which translates into a "thermal mass" limit 

Mt ~ l^^j M,. (24) 

The local condition for gap formation is presumably given by the larger of these two masses. 

Figure 7 shows radial profiles of My and Mt in three representative T-Tauri a-disk models 
for which we have estimated migration times previously. In each panel, the solid line corresponds 
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to the viscous mass, My, and the dashed Hne to the thermal mass, Mp. The upper panel shows 
the conditions for gap formation in our reference T-Tauri disk model (M = W~^Mq yr~^ and 
a = 0.02), while the middle and lower panels show these conditions in a high accretion rate model 
(M = W-'^Mq yr"^) and in a low a model [a =10 respectively. 

In each case, one easily identifies the strong dependence of My and Mr on the local disk 
aspect ratio, h/r, and how the location of the plateau at ~AU distances varies somewhat from 
model to model. For both the reference and high M models, the viscous mass. My, determines 
everywhere the conditions for gap formation, while the thermal mass, Mp, is relevant everywhere 
in the low a model. Gap formation conditions in models with intermediate values of a would be 
determined by both the viscous and thermal masses, depending on location in the disk. 

Except in the very innermost disk regions, effective gap opening conditions correspond 

to masses in excess of lOM® and reach values larger than 100 Earth masses typically beyond 
0.1-0.3 AU. In addition, it should be noted that the viscous mass threshold (Eq. [23]) was 
determined from a 2D formulation of spiral density wave theory. Since a more realistic treatment 
reveals that 3D Lindblad torques are significantly weaker than 2D ones, a derivation of My 
accounting for the finite thickness of proto-planetary disks should lead to even larger values of the 
viscous mass for gap formation. Therefore, provided one is interested in proto-planets with masses 
significantly smaller than the values shown in figure 7, estimates of local migration times ignoring 
disk feedback should be reasonably accurate, except for the effects of corotation torques, and 
except insofar as the collective action of many small planets may modify the disk. 



5. Discussion 

Our detailed calculations of Lindblad torques for low-mass proto-planets embedded in 
steady-state T-Tauri a-disks reveal strong sensitivities to radial properties of the background 
disk and its finite thickness. While we have found these results to be robust within the set of 
calculations we performed, several important limitations should also be emphasized at this point. 

We have ignored corotation torques for lack of an adequate parametrization of their effects. 
Although Lindblad torques are usually thought to dominate migration rates (e.g. KP93), the 
corotation resonances may be even more sensitive to local conditions and to feedback than the 
Lindblad ones, as explained in §3.3. Very likely, corotation torques would alter our quantitative 
results at order unity but would reinforce our main qualitative conclusions, namely, that the 
migration rate is strongly modulated by opacity transitions and other radially localized features 
of a disks. This will have to be tested in future work. 

It is unlikely that our steady-state constant-a disk models are accurate representations of 
T-Tauri disks in general. Indeed, the nature of angular momentum transport in proto-planetary 
disks is not well understood, even though it determines their structure to a large extent. Gammie 
(1996) observed that MHD turbulence driven by the magneto-rotational instability (Balbus & 
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Hawley 1991; 1998) should operate in the innermost regions (r ^ 0.1 AU) of T-Tauri disks 
under standard thermal ionization conditions, but probably not in their colder, outer regions. He 
proposed instead that MHD-turbulent accretion at large radii proceeds only in surface layers that 
are sufficiently ionized by cosmic rays. Glassgold et al. (1997) later refined this layered accretion 
picture, by taking into account the ionization due to X-rays from the central star. Since then, it 
has become clear that the radial extent of a layered accretion region is sensitive to details in the 
disk models (Promang et al. 2002; Matsumura & Pudritz 2003), with weakly-turbulent properties 
that were not previously anticipated (Fleming k, Stone 2003). 

Absent any contribution to angular momentum transport from hydrodynamical processes 
alone (see Balbus 2003 for a discussion), proto-planetary disks may also have to rely on planetary 
torques for their "viscous" evolution (Larson 1990). Goodman & Rafikov (2001) studied the 
propagation and dissipation of spiral density waves in gaseous disks, within the framework of 
two-dimensional theory. They showed that dissipation from wave breaking is quasi-local and 
they estimated a magnitude of angular momentum transport a ~ 10~^-10~^ for a uniform 
population of Earth-mass proto-planets embedded in a minimum mass solar nebula. It remains to 
be seen what efficiencies of angular momentum transport can be reached self-consistently in more 
realistic proto-planetary disk models that include explicitly the integrated torque feedback from a 
population of embedded proto-plancts. 

Clearly, these large uncertainties on the structure of proto-planetary disks will have important 
effects on migration time estimates. We believe, however, that the identification of opacity 
transitions in the disk as regions where migration rates vary rapidly with radius is not unique to 
our steady-state a-disk models. In the specific layered accretion model of Gammie (1996), for 
instance, opacity transitions cause variations in the slopes of radial solutions, much as in our 
models, and this may be a general property of models based on an a-formalism. 

What is less clear perhaps is the magnitude of this effect on radial profiles of migration times. 
In particular, our models have assumed a direct dependence of the opacity on the local values of 
density and temperature in the disk. In the presence of substantial radial mixing, for instance due 
to disk turbulence, these transitions could be made more gradual and the amplitude of the rapid 
radial variations of migration times we have found may be reduced accordingly. The issue of the 
turbulent nature of the disk is particularly important in view of recent results from fully turbulent 
MHD simulations for embedded proto-planets (Nelson &; Papaloizou 2003b; see also Laughlin, 
Steinacker &; Adams 2003). These simulations show how turbulent density perturbations can 
dominate over the specific spiral density pattern associated with an embedded proto-planet and 
lead to a largely random migration process. It remains to be seen whether long-term migration 
trends will emerge from such simulations and whether they will be sensitive to global radial disk 
properties such as those exhibited by our a-disk models. If not, it may be that our results on 
radial variations of migration times will be mostly relevant to the "dead zones" of layered accretion 
regions, where the turbulence is expected to be comparatively much weaker (Fleming k. Stone 
2003). 
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Not surprisingly, many of the uncertainties we have just highhghted provide good motivations 
for additional work on the subject. It would be interesting to include in our disk models a 
treatment of layered accretion similar to that described by Gammie (1997) and study proto-planet 
migration in that context. Another important extension of the present work would be to include 
the feedback on the disk due to planetary torques. This would allow us to study more accurately 
the conditions for gap formation, as well as the global efficiency of angular momentum transport 
in the disk associated with a population of embedded proto-planets. Ultimately, we may thus 
be able to study with some confidence the coupled evolution of a proto-planetary disk and its 
population of growing planets. 

6. Conclusion 

We have presented detailed calculations of migration times for low-mass proto-plancts 
embedded in steady-state T-Tauri a-disks, based on an approximate treatment of two- and 
three-dimensional Lindblad torques. We have emphasized the strong sensitivity of local migration 
rates to details in the background disk structure and we have argued that regions of opacity 
transitions may in general be characterized by sudden radial variations of the migration rates. 
Localized regions with the slowest migration rates may thus be preferred sites for gravitational 
interactions between proto-planets. In some of our disk models, we have found regions with 
migration times in excess of ~ 10^ yrs for a IQM^ proto-planet, and within the framework of the 
core accretion scenario, this could allow for significant mass buildup of a gaseous envelope before 
the proto-planet is accreted by the central star. 
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Fig. 1. — Profiles of various quantities in our reference model of a T-Tauri a-disk accreting steadily 
at a rate M = lO^^M© yr"^, with a viscosity parameter a = 0.02. (a) Profiles of surface density, S 
(solid line), Rosseland-mean optical thickness, r (short-dashed line), aspect ratio, h/r (long-dashed 
line) and Toomre Qt parameter (dotted line) are shown for both an irradiated and a non-irradiated 
disk. At large radii, the irradiated disk is thicker, less dense and less subject to self-gravity, (b) 
Profiles of midplane temperature, Tc (solid line), and internal effective temperature, T^s (dashed 
line; from viscous dissipation only), are shown for the irradiated and non-irradiated disk models. 
The irradiation temperature, Tirr, is shown as a dotted line. At large radii, the irradiated disk is 
hotter and it becomes vertically isothermal. 
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Fig. 2. — Normalized Lindblad torques for a disk with power law radial profiles of surface density, 
S, and temperature, T, and a fixed aspect ratio, h/r = 0.07. The two-dimensional formalism of 
spiral density wave theory is adopted here and the results of Ward (1997a) are recovered with high 
accuracy (compare his Fig. 3). (a) Inner, outer and net torques (dotted, dashed and solid lines, 
respectively) are shown for a disk with a fixed temperature profile {I = 1). (b) Same for a disk with 
a fixed surface density profile (A; = 0). 
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Fig. 3. — Same as Fig. 2, but this time adopting a modified formalism which accounts for three- 
dimensional effects in the torque calculation (see text for details). Note the clear reduction in all 
torque values and the qualitatively different dependence on the surface density profile. 
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Fig. 4. — Inward migration time of an Earth-mass proto-planet as a function its location in our 
reference T-Tauri a-disk model, with M = 10~^Mq yr~^ and a = 0.02. Solid and short-dashed 
lines compare the results obtained with the 3D and 2D formulations, respectively, for an irradiated 
disk. The dotted line shows, using the 3D formulation, the results obtained for a non-irradiated disk 
model (differences appear only at large radii). The long-dashed line shows the results obtained 
with a modified 3D formulation, in which the gas is assumed to be vertically extended over only 
half the local disk scale height. Migration times obtained with the 3D formulation are significantly 
longer than those obtained with the 2D one and they are more sensitive to opacity transitions in 
the disk. Regions of maxima of the migration times may be preferred sites for interactions between 
proto-planets. 
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Fig. 5. — Inward migration time of an Earth-mass proto-planet as a function its location in steady- 
state T-Tauri a-disks with M = 10~^Mq jr~^ and different values of the viscosity parameter, a. 
The 3D formulation of Lindblad torques is adopted and disk irradiation is included. Solid, dashed 
and dotted lines correspond to a = 2 x 10~^ (reference model), 5 x 10~^ and 10^'^, respectively. 
Decreasing the value of a increases the disk mass (which reduces migration times) and increases 
somewhat the values of radii at which opacity transitions occur. 
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Fig. 6. — Inward migration time of an Earth-mass proto-planet as a function its location in steady- 
state T-Tauri a-disks with a = 0.02 and different values of the mass accretion rate, M. The 
3D formulation of Lindblad torques is adopted and disk irradiation is included. Solid, dashed and 
dotted lines correspond to M = 10^^ (reference model), 10^'' and 2.5 x 10^*^ Mq yr^^, respectively. 
Increasing the value of M increases somewhat the disk mass (which reduces migration times) and 
increases the values of radii at which opacity transitions occur. 
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Fig. 7. — Radial profiles of viscous (solid line) and thermal (dotted line) masses above which a 
gap is formed (see text for details), in three of our steady-state T-Tauri a-disk models. The upper 
panel corresponds to the reference model (M = lO^^M© yr~^ and a = 0.02), the middle panel to a 
high-accretion rate model (M = lO^^M© yr"-*^ and a = 0.02) and the lower panel to a low a model 
(M = 10~^Mq yr"-*^ and a = 10~^). The viscous mass is generally the relevant one, except in low 
a disks. Both masses are strong functions of radius and they typically exceed 100 Earth masses at 
distances 0.1-0.3 AU. 



